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Abstract. Driven many-particle systems with nonlinear interactions are known to often display multi- 
stability, i.e. depending on the respective initial condition, there may be different outcomes. Here, we study 
this phenomenon for traffic models, some of which show stable and linearly unstable density regimes, but 
areas of metastability in between. In these areas, perturbations larger than a certain critical amplitude will 
cause a lasting breakdown of traffic, while smaller ones will fade away. While there are common methods to 
study linear instability, non-linear instability had to be studied numerically in the past. Here, we present 
an analytical study for the optimal velocity model with a stepwise specification of the optimal velocity 
function and a simple kind of perturbation. Despite various approximations, the analytical results are 
shown to reproduce numerical results very well. 



PACS. 89.40.Bb Land transportation 
83.60.Wc Flow instabilities 



45.70.Vn Granular models of complex systems; traffic flow 



1 Introduction 

While the field of traffic modeling and traffic simulation 
' has a long history [I], in the 90ies [2HI], it has also be- 
come interesting for a large community of physicists. Since 
then, various physical methods have been applied, rang- 
ing from fluid-dynamic and gas-kinetic approaches [M3 
[6] over many-particle models [10J up to cellular automata 
[9] . An overview of the respective literature has been given 
by various reviews [l0llll[[T2] . 

Besides computer-based studies, physicists have par- 
ticularly contributed with systematic (e.g. gas-kinetic) 
derivations [6] and analytical investigations [13] (see Ref. 
|10j for an overview). This includes the study of instabil- 
ties leading to a breakdown of free traffic flow [14,15,16, 
rT7in~8lfT9] , which has been studied experimentally only re- 
cently [21]. (Of course, empirical studies [2"2"]l2"3] have been 
already carried out for a longer time). 

While there are also models, where the instability of 
traffic flow depends on noise [2], we will focus here on 
models with a deterministic instability mechanism. Many 
traffic models become unstable in a certain density range, 
because of delays in the adaptation to changing traffic con- 
ditions. To avoid accidents, these delays are compensated 
for by over-reactions, which can increase a small initial 
perturbation and finally cause a breakdown of traffic flow. 
This kind of dynamics usually occurs in the medium den- 
sity range, where the change of the "desired", "optimal", 
or "equilibrium" velocity with a change in the distance or 
density is larger than a certain instability threshold [217] , 



To characterize the behavior of traffic flows more sys- 
tematically, Kerner and Konhauser have applied the no- 
tion of critical densities and critical perturbation ampli- 
tudes. Primarily based on numerical studies [24] . they 
found the following for a macroscopic, Navier-Stokes-like 
traffic model [10]: Altogether, there are four critical densi- 
ties p c k with k £ {1, 2, 3, 4}. Below some density p c \ > 0, 
any kind of disturbance eventually disappears. Between 
the densities p c ± and p C 2, one wide traffic jam builds up, 
given a large enough perturbation. A series of traffic jams 
may appear in a density range between p C 2 and some den- 
sity Pc3- Finally, a so-called "anticluster" can be triggered 
[24j , if the average density p is between p c % and some criti- 
cal density p C 4, while any disturbance disappears in stable 
traffic above p C 4. Similar observations have been made for 
other macroscopic traffic models [20J , but also microscopic 
ones of the car- following or cellular automata type [25, 26, 
[27] . The critical densities p c k depend mainly on the choice 
of the model parameters, in particular the relaxation time 
and the velocity-distance or velocity-density relation. 

If the average density p falls into the density ranges 
\pciiPc2] or [pc3,Pvi]) traffic flow is predicted to be 
metastable, i.e. characterized by a critical amplitude 
Ap cr (~p) for the formation of traffic jams. This amplitude 
is zero for p = p C 2 and p = p C 3, i.e. at the boundaries 
of the linearly unstable regime, while the critical ampli- 
tude grows towards the boundaries p c \ and p C 4 of sta- 
ble traffic and does not exist beyond these values. Per- 
turbations with subcritical amplitudes Ap < Ap CT (jj) are 
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Fig. 1. (Color Online) Formation of a traffic jam in case of 
overcritical perturbations (blue lines) and relaxation to free 
traffic flow in case of undercritical perturbations (red lines) 
for a low vehicle density. The plot shows 10 trajectories rep- 
resenting each 10th vehicle. Altogether, the simulation was 
performed with N=100 vehicles on a circular road of length 
L — N/p. Time is measured in units of the relaxation time r, 
vehicle locations in units of the safe distance do- The slope of 
the lines corresponds to the vehicle speeds in units of cLq/t, and 
their vertical separation reflects the vehicle distances in units 
of do- Therefore, the blue lines correspond to stop-and-go traf- 
fic, where the vehicle distances in the stop regime are small, 
i.e. vehicle density is large. The red trajectories show that the 
same number of vehicles can constantly move at the free speed 
v°, if the initial perturbation is small enough. 



eventually damped out (analogous to the stable density 
ranges) , while perturbations with supercritical amplitudes 
Ap > Ap cl (jj) grow and form traffic jams (similar to 
the linearly unstable density ranges). The situation in 
metastable traffic is, therefore, similar to supersaturated 
vapor [24] . where an overcritical nucleus is required for 
condensation ("nucleation effect"). 



Fig. 2. (Color Online) Formation of stop-and-go traffic from an 
initial perturbation in linearly unstable traffic flow. For details 
see the main text and the caption of Fig. Q] 



In order to gain analytical insights into the mecha- 
nisms of nonlinear instability, we will proceed as follows: 
In Sec. [2 we will introduce the optimal velocity model we 
will work with. Afterwards, in Sec. [3j we will discuss how 
perturbation sizes can be measured and what are neces- 
sary preconditions for the existence of traffic jams. We will 
then determine sufficient conditions by calculating critical 
densities and critical perturbation amplitudes for a simple 
kind of perturbations: In Sec. HJ all vehicles on a circular 
road but one are assumed to be in a big traffic jam, while in 
Sec. [3J all vehicles but one are assumed to experience free 
flow conditions. Finally, Sec. [6] will summarize our results 
and present an outlook. Supplementary, and for compari- 
son with our results, the Appendix gives a short derivation 
of the characteristic constants of the car-following model 
studied in this paper. 



Fig. 3. (Color Online) Formation of a traffic jam in case of an 
overcritical perturbation (blue lines) and drop of the vehicle 
flow to zero in case of undercritical perturbations (blue lines) 
for a large vehicle density. For details see the main text and 
the caption of Fig. [l] 



2 Introduction of the Applied Car-Following 
Model 

In the past, a large number of papers has addressed the 
instability of traffic models in an analytical way, and our 
discussion naturally needs to restrict itself to a few of 
them. For example, Kerner et al. [28] have presented an 
asymptotic theory of traffic jams for a Navier-Stokes-like, 
macrocopic traffic model, while instability analyses for 
car-following models were carried out by Herman et al. 
[14] , Bando et al. |15j . and others (see citations in Ref. 
[TU]). Some recent analyses of traffic instabilities are by 
Mitarai and Nakanishi [T6] and by Wilson [1711181119] . The 
most relevant references on which this paper builds on are 
Refs. [2"9"1[2"T] . For a car- following model and a related cel- 
lular automaton, respectively, these have presented ana- 
lytical calculations of the characteristic constants of traffic 
flows such as the dissolution speed c of traffic congestion 
(which agrees with the propagation speed of traffic jams), 
the outflow Q ou t from congested traffic, or the vehicle den- 
sity Pja.ni in congested traffic for a car-following model and 
a related cellular automaton, respectively. 

The success in determining the self-organized con- 
stants of traffic flow analytically is based on boiling down 
vehicle dynamics to its essence. That is, we will work with 
a very simple model (a "toy model"), that is not aiming 
to be fully realistic, but instead at performing analyti- 
cal calculations and deriving stylized facts that would not 
be analytically accessible for a more realistic model with 
a reasonable amount of effort. The model we will focus 
on is the car-following model by Bando et al. [15] with 
a stepwise optimal velocity function. It is defined by the 
equation of motion dxj/dt = Vj(t) relating the change of 
location Xj (t) of a vehicle j in time t with its speed Vj (t) 
and the acceleration equation 



dvj v (dj(t)) - Vj (t) 



dt 



(1) 



The parameter r > corresponds to a small "relaxation 
time" . The "optimal velocity function" v (d) depends on 
the vehicle distance dj (t) = 
its predecessor ( "leader" ) j 



Xj-i{t) — Xj(t) of a vehicle to 
- 1, and it is specified here as 



v a (d) 



v° if d > d Q , 
otherwise. 



(2) 



do > v°t is a safe vehicle distance and v° the free speed 
(or speed limit) of the vehicles, which are assumed to be 
identical, here. Defining the local vehicle density p as the 
inverse of the vehicle distance and the stationary and ho- 
mogeneous equilibrium flow as Q c (p) = pv (l/p), we find 



Qc( P ) 



pv° if p< l/d , 
otherwise. 



(3) 



It is obvious that this flow-density relationship ("funda- 
mental diagram") is not realistic. However, when the traf- 
fic flow is unstable with respect to perturbations in the 
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flow, a much more realistic, self-organized flow-density re- 
lation results, namely the so-called "jam line" [28] 



J(P) 



1 



Pjam-' 



(Pjam - P) 



(4) 



(see Fig. 3]). The jam density pj am and the delay T between 
the acceleration of successive vehicles at the jam front can 
even be calculated (see Ref. [29] and the Appendix of this 
paper). The corresponding implicit relations are 



ft 



— = d Q - v°t(1 - e~ T / T ) 



(5) 



and 



T = 2r(l-e- T / r ) . (6) 

It is also worth stating that the characteristic outflow 

1 



Qout — PoutV 



(7) 



d /v° + T/2 

from congested traffic is smaller than the maximum flow 

„0 



do ' 



(8) 



that is, there must be a density region in which two differ- 
ent flows are possible. Hence, the actually assumed density 
value depends on the history. This is called a hysteresis ef- 
fect. Once a traffic jam forms, there is an effective capacity 
drop of size 



^^drop ^;n 



1 



V°T/2 



d (d /v° + T/2) 



>0 



(9) 



(see Fig. E]). 



Fig. 4. Schematic illustration of the traffic flow as a function of 
the vehicle density p. The density ranges between p c i and p C 2 
and between p c z and p C 4 correspond to metastable traffic flow 
(see main text). The quantities Q c k = Q(Pck) denote the flows 
belonging to the critical densities p c k, which will be analyti- 
cally calculated in this paper. The thin solid line corresponds 
to the jam line J(p). Its intersection point with the free branch 
of the fundamental diagram defines the value of the character- 
istic outflow Qout and the density p ou t, while its intersection 
point with Q = defines the jam density Pj am . The fundamen- 
tal diagram, i.e. the flow-density relation Q c (p) in the station- 
ary and homogeneous case, is represented by thick solid lines. 
The difference between the maximum flow Q ma x = max p Q c (p) 
and the characteristic outflow Qout is called the capacity drop 

<!-ittf drop ■ 



3 Definition of Critical Amplitudes and 
Critical Densities 

It is very encouraging that characteristic constants of 
traffic flows such as the jam density Pjam, the outflow 
Q out from congested traffic, and the propagation speed 
c = — l/(pj am T) of jam fronts can be analytically derived 
from the optimal velocity model. This just requires a sim- 
ple enough specification |(2]) of the function v (d), see Ref. 
i29j and the Appendix of this paper. Therefore, we will 
try in the following to derive the critical densities and 
critical amplitudes for this model as well. For this, let us 
study a circular one-lane system of length L with N ve- 
hicles and an average vehicle density p = N/L. In order 
to avoid finite size effects, we will assume a large system 
with many vehicles. According to Eqs. ([T]) and ^j, if all 
vehicles have a distance greater than do to their leader, all 
of them will accelerate and reach the maximum speed v° , 
while we expect congested, standing traffic, if all vehicles 
have a distance smaller than do, since this forces them to 
decelerate. 

Therefore, an interesting dynamics will only occur if 
some vehicles have distances larger than do, while others 
have distances smaller than do- In the following, we will 
focus on this case. One first observation is that linearly un- 
stable traffic should be possible only for p = N/L = 1/do, 
so that, for the stepwise specification |2]) of the optimal 
velocity function v a (d), we have the special case 



1 



Pc2 = Pc3 = 



(10) 



But what is the value of the critical density p c \ at which 
non-linear instability starts to be possible, and is there a 
critical density p C 4, beyond which even arbitrarily large 
perturbations will fade away? And can we approximately 
determine the critical amplitudes Ap cr (j))l 

Let us in the following focus on a special kind of per- 
turbations (see Figs. [5]and [7]): We will assume that (JV— 1) 
vehicles have an identical distance d/ ) to their respective 
predecessor, while the distance di(0) of the first vehicle to 
its predecessor is 

di(0)=£-(i\r-l)d (O ). (11) 

Accordingly, we define the perturbation size as 



Ap = 



1 



1 

^^rfiTo) 

L-Nd 



1 
i(0) 



1 



L-(N- l)d ( o) 



(0) 



d {0) [L-(N-l)d i0) ] 



(12) 



In order to avoid finite-size effects, we will assume a very 
large system of length L — N/~p with N S> 1 vehicles. 
Then, we have 



d(o) = 



L-di(0) 



JV-1 



L 

N 



which implies 



Ap: 



di(0) 



(13) 
(14) 
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The critical amplitudes Ap CT are basically defined by not 
changing in time (neither growing nor shrinking), i.e. by 
marginal stability. Moreover, the critical density p c \ is 
characterized by the fact that the critical amplitude ceases 
to exist for lower densities, and the same applies for den- 
sities greater than the critical density p C 4. 

Remember that the propagation of perturbations re- 
quires at least one vehicle distance (either d(n) ~ 1/p or 
di(0)) to be above and another one below do- This allows 
us to estimate a lower bound for the critical perturbation 
amplitudes Ap CI (~p). Two cases my be distinguished: 

1. If the density ~p ss l/<i(o) is greater than 1/do, we must 
havedi(O) > d , which implies [p— l/di(0)] > p— 1/do- 

2. If the density p « l/^(o) is smaller than I /do, we must 
havedi(O) < d , which implies [l/di(0) — p] > 1/do— ~p. 

Consequently, the critical amplitude can become zero only 
for p = 1/do = Pc2 = Pc3, and altogether we must have 



A°cr(p) > 



1 

do 



(15) 



However, this is only a necessary condition, while we will 
determine sufficient conditions for the existence of stop- 
and-go waves in the following sections. In our further anal- 
ysis, we will treat the case of large densities separately 
from the case of small densities. 



4 The Case of Large Densities 

The situation for densities p > 1/do is illustrated in Fig. 
[5j We assume that all vehicles j start at time t = 
with speed 1^(0) = 0. Moreover, all vehicles but one are 
assumed to have the initial distance di ) < do, while 
the remaining vehicle j = 1 has the distance di(0) = 
L — (N — l)d(o) > d , where L is the length of the as- 
sumed circular road. It will turn out that the simplicity of 
this initial perturbation is the reason for the feasibility of 
our calculations. Without loss of generality, vehicle j — 1 
starts to accelerate at x = and t = 0. 

We will now have to identify the possible reason for 
a disappearance of the initial perturbation in the course 
of time. In the case of large densities p > p c ^ = 1/do, an 
initial perturbation will fade away if the maximum vehicle 
distance does not allow to reach a sufficiently high speed in 
the acceleration process to reach a vehicle distance equal 
to or smaller than d( ) after its successive braking maneu- 
ver. While the speed in the vehicle queue is approximately 
zero (despite for, maybe, the last few vehicles in the queue 
which may still decelerate), the maximum speed V(o) is 
reached at the time trn\ when the vehicle under consider- 



<"(0) 

ation starts to decelerate. Focussing on vehicle j 
find 

v i o)=vi(t { o ) )=v°(l-e' t ^/ T ), 



1, we 



(16) 



Fig. 5. Illustration of the initial distribution of vehicles on a 
ring road assumed in the high-density case (p > 1/do). Vehicles 
j € {2, . . . , iV} are standing in a traffic jam with speed vj — 
and a distance d(p) < do to the predecessor. Vehicle j = 1 starts 
with zero velocity as well, but can accelerate, if its distance 
di(0) to the next vehicle is greater than do. 



As the first vehicle's position at time t = is assumed 
to be xi(0) = 0, the time-dependent location of vehicle 1 
is given by 



Xl (t) = / dt'vi{t') = / dt' (v° - w°e-*' /r ) 



= w t + w°r(e- t / T -l). (17) 

The time t( ) at which vehicle j — 1 starts to decelerate is 
the time t at which its distance to vehicle j — N, i.e. the 
last vehicle in the queue, becomes do- Since vehicle j = N 
is located at xn = L — (N — l)d(o), this implies 

L-(N-l)d i0 )-xi(t {0) )=do (18) 

or 

L-(N- l)d (0) - d = v°(t {0) -t + re-W T ) , (19) 

which is an implicit equation determining the acceleration 
time period tr )- At the end of its deceleration process, the 
previously first vehicle of the queue, which has then joined 
the end of the queue, will have a distance d(i), which de- 
pends on the maximum speed U( ) • Since the deceleration 
process according to Eqs. (p} and ([2]) obeys an exponential 
velocity decay 

v l (t)=v (0) e-( t - t w)/T, (20) 

the resulting minimum distance can be determined as 

00 
d ( i) -d - J dt z^e-C-W/- = do - v {0) t . (21) 

*(0) 

The other vehicles in the queue are expected to have the 
same distance to their respective predecessor after one cy- 
cle of acceleration and deceleration. That is, while their 
distance was dm\ in the beginning, their distance after one 
cylce will be dm , after two cycles it will be d^) > and so on. 
The iterative equations to determine the decisive quanti- 
ties can be derived analogously to Eqs. ([2Tj) and (fl6|) : 



d (n+1) = d - v (n) T = do - v°t(1 - e-*<")/ r ) , (22) 



where ti n \ as a function of di n \ is determined from 



because Eqs. ((T|) and @ imply v\(t) = v°(l— e t ^ T ), which 
can be easily checked by differentiation with respect to t. 



L-{N-l)d {n) -do = v a (t {n) 



T + TC 



-t(„)/T 



(23) 
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similarly to Eq. (fT9| . Equation (|23|) allows us to replace 
the right-hand side of Eq. ([22]) , which yields an equation 
for d( n +i) as a function of dt n y. 



d( n +l) = L- (N - l)d(„) - V°t(„)(d( n )) . 



(24) 



It is, therefore, interesting to ask, whether the series d( n ) 
converges and, if yes, to what value. If the values d/ n \ stay 
the same for different values of n, the initial perturbation 
is stable over time. If the distances go down, then the 
perturbation grows. However, if the values of di n ) grow 
with n, the initial perturbation fades away (which is ex- 
pected to happen, when the perturbation is too small). 
Therefore, the critical amplitude is given by the condition 
d( n +i) = dr n \ of marginal stability, which together with 
Eq. d24j) implies d( n ) = L - (N - l)d(„) - u°t( n ) or 



V°t(n) 



L- Nd 



(n) 



(25) 



Multiplying Eq. (25]) with (TV - 1), we obtain for the 
marginally stable case d( n +i) = d( n ): 



(N - l)d (n) = (N-1) [d - «°r(l - e-*(-)/ r ) 
Inserting Eq. (|23|) finally yields 

u%„) = i - ATd + AW(l - e~'<"> /r ) 



(26) 



(27) 



Dividing Eq. (J2T[) by AT, considering L/iV = 1/p, and per- 
forming the limit N — > oo gives 



^ /T )=do-~, 
P 



«°r(l- 

i.e. one solution with a finite value of t(n) 
do - 1/p 



H") 



= -rln 1 - 



— ( do - = 
v" \ p 



(28) 



(29) 



where we have applied d( ) ~ 1//3 according to Eq. (TH 

Equation (|32p agrees well with our numerical findings 
(see Fig. [6]). Obviously, the critical perturbation ampli- 
tude is zero for p = 1/do = Pc3, as it should. Moreover, 
for the specification (TJJ) of the optimal velocity function, 
the critical amplitude does not diverge at a finite density. 
Nevertheless, there is a critical density p c ^, which is given 
by the fact that d(„ + i) > dp — v°t according to Eq. 



Considering di n 



+i) = d(o) ~ 1/p in the marginally sta- 
ble case, a critical amplitude ceases to exist for densities 
larger than 

P ci = 1 b~ • ( 33 ) 

do — v U T 



Fig. 6. (Color Online) Results of computer simulations of the 
car-following model Jl} with the optimal velocity function ([2} 
for a time discretization St = 0.1 and v° = do/r (where do = 1 
and t — 1 has been assumed because of the possibility to scale 
space and time). The simulations were run for N — 100 ve- 
hicles, with the initial conditions illustrated in Figs. [5] and [7] 
Given a certain average density p, the length L of the simu- 
lation stretch was chosen as L = N/~p, and the perturbation 
amplitude Ap CT (p) — |l/d(o) — l/[£ — (N — l)d(o)]| determined 
the initial distance d(p) between TV — 1 of the vehicles. The bold 
dashed lines represent the necessary condition (|15|) for the exis- 
tence of stop-and-go waves, while the thin solid lines represent 
the approximate critical amplitudes calculated in this paper. 
The critical amplitudes according to our numerical simulations 
lie between the areas represented by different kinds of symbols 
reflecting different traffic states at t — 2000: green plus signs 
mean all vehicles move at the desired speed v°, blue circles 
mean all vehicles have stopped moving, and red dots mean ve- 
hicles have different speed values, corresponding stop-and-go 



Here, our approximation is based on the first-order Taylor 
expansion ln(l— x) « —X. Together with Eq. ([25]) . we have 



do- - 
P 



v°t 



(n) = L - Nd( n ) 



(30) 



and considering d( n+1 ) = d(„) = . . . = dm) finally gives 



L- Nd 



(o) 



do - r 

P 



(31) 



This can now be inserted into Eq. (TT2"]) to obtain the crit- 
ical amplitudes: 



^Pcr(p) 



L- Nd 



(0) 



d (0) [i-(iV-l)d ( o)] 
do - 1/p 



d(o)(do - l//9 + d( )) 
do - 1/p 



do/p 



do 



(32) 



5 The Case of Low Densities 

The initial condition assumed for low densities is illus- 
trated in Fig. [7] While at large densities, vehicle j = 1 
had a larger distance than the other vehicles, it has now a 
smaller distance. In the following, we will focus on vehicle 
j = N, of which we assume that it is located at x = 
and starts to accelerate at time t = with an initial speed 
^iv(O) = 0, while the following vehicle j = 1 is assumed to 

have the initial distance di(0) = d m ( n < d . As the other 
vehicles have, by definition, the same initial distance d(o) 
from each other, the distances are related via the equation 



,(0) 



di(0) = d^ = L-(AT-l)d (O ) 



(34) 



Equation (|12|) determines again the perturbation size. In 
the case of small densities p < p C 2 = 1/do, the survival of a 
perturbation requires that there is at least one vehicle with 
a distance d min < do for a long enough time period t* to 
force the successive vehicle to brake. (If this vehicle would 
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already start to accelerate before the follower reaches a 
distance do to it, this would cause the perturbation to fade 
away.) Therefore, let us determine i* in the following. 



Fig. 7. Illustration of the initial distribution of vehicles as- 
sumed in the low-density case (p < 1/do). Again, all vehicles 
are standing at time t — 0, but vehicle j = 1 has the initial 
distance di(0) = a^[ n — L—(N—l)d(o), while all other vehicles 
have the initial distance dio) > do. The initial vehicle speeds 
were set to zero. 



5.1 Derivation Focused on the Trajectory of One 
Vehicle 

According to Eqs. (dJ and ^, the speed of vehicle j — N 
evolves in time according to dvjy/dt = [v° — i>jv(£)]/t, 
which implies 



< x ! / 1 = / dt' ^jP- =v°(l- e"*/ T ) « «° , (35) 



and the related distance moved is 

t 
d N (t) = / dt' v{t') = v°t + v°r(e-* /T - 1) « v° ■ (t - r) , 

o 

(36) 

where the approximate equalities hold for f >t. Let aLL 
denote the distance of the vehicles j G {2, . . . , N—l}, after 
they were stopped by vehicle j — 1 or a follower. Then, 
the distance moved by vehicle j = N before it starts to 
decelerate at time t = t* (when it has reached a distance 
do to its predecessor), is 

d N (U) = L- d { ± (N 2)d^ n - do . (37) 

d lnin denotes the minimum distance of a vehicle at the end 

(2) 

of its first deceleration maneuver, d [ after the second 



for the survival of a perturbation. As the first deceleration 
maneuver extends over a time period At = t* — T(o) — 
(N - 2)T (1) , the evolution v N (t) = v e-^~ t ' 1) '>/ T of the 
vehicle speed for £; < t < ti + At implies that the 
minimum distance d„L afterwards is 



2 (1) 



mm 
t^ + At 

i - dt' v N (t') 



d - v°t(1 - e~ M/T ) >d -v°i 



(40) 



Moreover, let us assume a subsequent increase of the speed 
according to u°(l — e _ (*~* )/ r ), where t' represents the 
starting time of the acceleration maneuver. Then, since 
the following vehicle starts to accelerate when that vehicle 
has reached a distance do, the time period Tm between 
successive acceleration maneuvers is determined by the 
equation 



t'+T (1) 



do = d^ n + I *Al-«" M/T ) 



T, 



(i) 



r(l 



-T m /r 



(41) 



Consider now that, for the critical amplitude, the condi- 



tions d, 



i) 



d ( °L and T (1) 



must hold. Then, Eqs. ([39]) and (41]) imply 



T( ) of marginal stability 



(N - 1)T (0) = (N - 1) 



'a a w 

«0 - a min 



•r(l 



_ P - r to)A 



L-(N-l)d ( Z~do 



(42) 



Rearranging this and dividing the result by N, in the limit 
N > 1 and with L/N = 1/pwe get 



i - d w w°r(l 
P 



" T (0)/^) 



(43) 



Furthermore, with d^ n = d^ n and T (1) = T (0) , Eq. gTJ 
leads to 



' nun //T-i \2 

one, etc. Equations ((36]) and ((37]) imply the start of the ^(°) _ j _ v " T n - e _T < > /T ) -v°T (0) « u° (0) , 



deceleration maneuver at time £ = £* with 



t (i) =T+ L - d£j - (N - 2)d^ n - dp _ g 



In order for the perburbation to persist, this time period 
must be shorter than the time period Tr \ + (N — 2)Tm, 
at which vehicle j = 1 starts to accelerate, where Ti n \ 
is the time shift between the acceleration of two subse- 
quent vehicles standing at a distance d^" n . This implies 
the threshold condition 



(44) 



where the last approximation is based on a Taylor expan- 
sion of second order. Therefore, we have 



(«°T (0) ) 2 a 2™° (d 



7(0) 



(45) 



(see Fig. H]), while Eqs. (O and ((32]) together imply 



v°T {0) * - - d[, 
P 



(0) 



(46) 



r + L (N l )dmin ^ = r ( o) + {*- 2)r ( i) (39) 



As a consequence of Eqs. ([4"5"j) and (05]), we have 

1 ) 2 «2r« (d -dil). (47) 



- - d (0) 

— "min 
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Solving this with respect to d m ( n finally gives 



,(o) 



± 



/(™S 



2TV°d 



but only the solution with the plus sign meets the require- 
ment d m ( n > in the relevant density range (see Fig. [9|) . 
Considering Eqs. (fT2|) and (fl3|) . we find the following re- 
lationship for the critical amplitude: 



Ap c 



1 

ho) 



l 



di(0) 



1 



d (0) 

"min 



(49) 



(0) 

min 



The critical amplitude ceases to be well-defined, when d, 
assumes complex values, i.e. when the expression in Eq 
(j4"8|) under the root becomes negative. The critical density 
corresponds p c \ to the density for which the root becomes 
zero, i.e. both solutions in Eq. (j48| agree. This leads to 



1 



Pel = 



d Q + tv°/2 



(50) 



With this, we have successfully estimated the critical den- 
sity and the critical amplitudes, as Fig. [6] shows. Notably, 
our results for the low-density regime deviate significantly 
from the lower limit given in Eq. (fl5|) . It is for the first 
time that the problem of determining critical perturbation 
amplitudes of a traffic model has been analytically solved. 



Fig. 8. (Color Online) Exact solution of Eq. (|5ip as compared 
to its approximate solution (|45[) . The approximation works par- 
ticularly well in the relevant range of d\J ixx /do ~ 1. 



Fig. 9. Both real- valued solutions of Eq. (|48[) . It can be seen 
that the lower solution would predict negative values of d m ( n for 
most density values p. Therefore, the upper solution must be 
chosen. Below p = Pel, the solution of Eq. (148 [I has an imag- 
inary part, which does not have a reasonable interpretation, 
here. 



5.2 Derivation Focused on A Leader and Its Follower 

While we have based our previous calculations on the 
study of vehicle j — N, we may also derive our results 
from the consideration of two successive vehicles instead. 
Analogously to Eq. (J4I 
j = 1 starts to accelerate is given by 



the time T( ) after which vehicle 



7 (0) 



(0) 



Tl 



a - T to)/ 



T ) 



d - 



(51) 



At that time, the following vehicle has the distance 



d ( X = d i0 )-v°t" 



-v°t(1 



<T {0) -t")/7 



where 



d(o) - dp \jp - dp 



(52) 



(53) 



is the time, at which vehicle j = 2 starts to decelerate, 
while it moves at speed v° before. For the survival of the 
perturbation we have to demand t" > and d\ I J n < d„ lin . 
In first-order Taylor approximation, we have 



d (] 



7(D 



v°t(1 



e -(T( )-t")/r 
*V°(T {0) -t") 

so that together with Eq. (f5"2"j) we get 



v ° T (o) ~ I = - <k 



v"T< 



1 



(0) = - 
P 



7(D 



(54) 



(55) 



As pointed out before, for the critical amplitude we can 
presuppose marginal stability with cr_j_ = di?. Consc- 



7(0) 



results 



quently, the relation (w T( )) 2 ~ 2rv°(dp 

again as an approximate solution of (|?T|) . see Fig.fTUl 



Fig. 10. Critical amplitude according to Eqs. (T35J and (fiSJ). 
The curve is very well compatible with the numerical results, 
as shown in Fig. [6] 



6 Summary and Conclusions 

Our findings can be summarized as follows: Based on the 
car-following model ([T]) with the stepwise optimal velocity 
function ^ and a simple type of initial perturbation, it 
is possible to derive critical amplitudes and critical densi- 
ties. Linearly unstable behavior, where the breakdown of 
free traffic flow is caused by minor perturbations, is found 
only between the critical densities p c i and p c ^. In case of 
the stepwise optimal velocity function ((2]), both critical 
densities agree, and we have p C 2 = p C 3 = 1/cfo. Below this 
density, the survival of a perturbation needs at least one 
vehicle with a distance d m [ n < dp for a long enough time 
period to force the successive vehicle to brake. This re- 
quires densities between the critical densities p c i and p C 2, 
where p c i is given by Eq. ([50)1 in good agreement with nu- 
merical results, see Fig. [6) Furthermore, the perturbations 
must have a size greater or equal to the critical perturba- 
tion amplitude Ap CI (p). Equation (|4"9"| with (|4"g)l represent 
the related results of our approximate analytical calcula- 
tions. 

If the average density p exceeds the value 1/do, an ini- 
tial perturbation will fade away, if the maximum vehicle 
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distance does not allow drivers to reach a sufficiently high 
speed in the acceleration process to reach a vehicle dis- 
tance equal to or smaller than d( ) after the successive 
braking maneuver. With Eq. ([33|) . we have found an an- 
alytical formula for the critical density p c i, above which 
perturbations will necessarily decay. However, for average 
densities between p c s and p C 4, perturbations larger than 
the critical amplitude given by Eq. (|32[) will grow and form 
moving jams of the kind of "anticlusters" , while smaller 
perturbations will fade away, giving rise to homegonous, 
congested traffic (and in the special case of the stepwise 
optimal velocity function assumed here, even to standing 
traffic). Note that, for other specifications of the optimal 
velocity function, homogeneous congested traffic flows will 
usually be finite. Furthermore, it should not be forgotten 
that the critical perturbation amplitude may depend on 
the shape of the initial perburbation. 

Our analytical findings are summarized by Fig. [6l 
Note that the self-organized flow-density relation (the 
"jam line") differs significantly from the "fundamental di- 
agram" resulting from the optimal velocity function in the 
case of stationary and homogeneous traffic flows, which, 
however, may be unstable. An interesting observation for 
the stepwise specification of the optimal velocity function 
is that large perturbations at high vehicle densities can 
reach greater average flows than small perturbations, while 
the situation at low densities is characterized by a capacity 
drop, see Eq. ©. 
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calculated (see Ref. [29]). In the following, we summarize 
the calculations in our notation, here. For this, let us as- 
sume a number of jammed vehicles with velocity zero and 
distance smaller than do. Pjam shall denote the jam den- 
sity. 

If we assume that a car starts to accelerate out of the 
jam only when its leading car has approximately reached 
its desired velocity v°, and starts to decelerate when its 
predecessor has almost stopped in the traffic jam, we have 
the acceleration equation dvj/dt « [v° — Vj(t)]/r > for 
dj > d and dvj/dt « [0 — Vj(t)]/r < for d < d . Conse- 
quently, we find 



Vj(t) 



>(l_ e -(*-*o)/ T ) if dj{t) >d , 



V e -(t-t 2 )/7 



otherwise. 



(56) 



to is the time point when the acceleration of vehicle j 
starts with with dj (to) = do and ddj (to)/dt > 0, while £2 is 
the successive time point with ^(£2) = do and ddj(t2)/dt < 
0, when the deceleration starts. 



Fig. 11. Illustration of the approximate periodic change of dis- 
tance and speed ("hysteresis loop") for fully developed traffic 
jams (solid lines with arrows), after Ref. [29]. The dashed line 
indicates the underlying optimal velocity function v (d) , which 
is very different. 



This implies a linear increase of distance with speed 
(see Fig. ITTj) . Consequently, the desired velocity 
Vj(t) = v° is reached at the distance 



=do+TV°(l 

Pout 



-T/t\ 



(61) 



which defines the density p out related with the charac- 
teristic outflow Q ou t = Poutv from congested traffic. 
Phase 2 is characterized by dVj-i(t)/dt < and 
dvj(t)/dt > 0, i.e. vehicle j — I already decelerates, 
while vehicle j still accelerates. This phase is assumed 
to start at time t\ > to, and we have ddj(t)/dt = 
yOg-Ct-tO/r _ v 0(i _ e -(*-t )/T) < Therefore, we 
get 

dj(t) = dj(h) - u°r(e- ( *- tl)/T - 1) - v° ■ (t - h) 



v°T(e- 



-(t-t )/r 



-(ti-to)M 



« dj(h) - v°(t - h) + v°r(l - e-t*-*^ , (62) 

where the approximation assumes t — to > ti — (o>t. 
That is, the distance dj(t) goes down, while the speed 
Vj(t) — v°(l — e _ ' <_t °''' T ) is approximately v°. The 
distance dj(t) becomes do at time £2 with 



(*2-*l) 



dj(h) -d 



(63) 



We may distinguish four different phases: 

Phase I is characterized by dvj-i(t)/dt > and 
dvj(t)/dt > 0, i.e. vehicle j and its leader j — I both 
accelerate. According to the travelling wave concept, 
the trajectory of vehicle j is exactly identical with the 
one of its predecessor, but shifted by some time period 
T, which corresponds to the delay between the accel- 
eration of two successive vehicles out of the traffic jam, 
i.e. 



«*(*) 



Vj-l 



(t-T). 



(57) 



The speed c of jam resolution is then given by the dis- 
tance — 1/pjam between two jammed vehicles, divided 
by this time period T: 



(58) 



Pjam-* 



The negative sign is a consequence of the fact that 
vehicle j stands upstream of vehicle j — 1, but accel- 
erates later, so that the downstream jam front travels 
opposite to the direction of motion. 
Considering the delay T, we get 



ddjjt) 

dt 



= tt,_i(£) -Vj(t) 
= Vj (t-T)- Vj (t)., 



(59) 



which by integration over time and with dj(to) = do 
results in 



dj(t) 



1-e 



>°r(f 
= d + (l-c- T / T )r^(£) 



-(t-to)A 



(60) 



Phase 3 starts at time £2 > £1 and is characterized by 
dvj-i(t)/dt < and dvj(t)/dt < 0, i.e. vehicle j — 1 
decelerates, and vehicle j does the same with a time 
delay of T. One can say that phase 3 is the inverse 
process of phase 1, and we get 



d 3 (t) = do -«°(1 - e- T / T )r{l - e -('-^)/-) 



(1 - c- t ")t[v -v s {t)] 



(64) 



Accordingly, the distance is monotonously decreasing 
with time. The minimum distance is reached at time 
£3 for the jam density Pjam, which is defined by 



1 



= do-TV°(l-C- T / T ) 



Pjam 

Together with Eq. ([BT|) we find 
I I 

Pout Pjam 



2tv°(1 



-T/t) 



(65) 



(66) 



Phase 4 starts at time £3 and is characterized by 
dvj-i(t)/dt > and dvj(t)/dt < 0, i.e. vehicle j still 
decelerates, while its leader j — 1 already accelerates. 
Phase 4 is the inverse process of phase 2, and we have 
ddj/dt = v°(l - e^-* 3 )/ 7 ") - v°e- ( - t - t2 ^ T . Therefore, 
we find 

dj (£) = d,(£ 3 ) + v° ■ (t - £3) + v Q T(e-^-^/ T - 1) 

+ U °r(e- (t - t2)/T -e- ( * 3 - t2)/T ) 

w — + v°(t - £ 3 ) - v°t( I - e -(*-* 3 )/ TN ) , (67) 

Pjam V ' 
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where the approximation assumes t — £3 > £3 — 1% 3> r. 
Therefore, the distance grows in time, while the veloc- 
ity Vj(t) — ti°e~(*~ t2 )/ r of vehicle j is approximately 
zero. At time £4, the distance dj(t) becomes do- The 
difference £4 — £3 determines the time delay T between 
two successive acceleration maneuvers of cars leaving 
the traffic jam. Inserting the definition (|65jl for /Oj a m, 
we obtain 



d a = — +v°T-Tv°(l-e- T / T ) 

Pjam V ' 

= d + v°T - 2tv° (l - c- T / T ) (68) 



and the implicit relationship for the time shift 

T = 2r(l-c- T/r ). (69) 



Together with Eq. (fBTj) . this yields 



1 


, v°T 





= d + —- 


Pout 


2 



(70) 



Therefore, according to Eqs. ([63]) . (|6"Tj) , and (f6"9"|) , it is 
possible to express the characteristic constants T, c, Pj am , 
and Pout through the model parameters do, v , and r only, 
without any dependence on the initial condition. Note that 
all the characteristic parameters are no model parameters. 
They are rather a result of the self-organization of char- 
acteristic jam fronts (see Fig. Hi]) . 



